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Purpose: Corneal epithelium is maintained by a population of stem cells (SCs) that have not been identified by specific 
molecular markers. The objective of this study was to find new putative markers for these SCs and to identify associated 
molecular pathways. 

Methods: Real time PCR (rt-PCR) was performed in 24 human limbal and central corneal epithelial samples to evaluate 
the gene expression profile of known corneal epithelial SC-associated markers. A pool of those samples was further 
analyzed by a rt-PCR array (RT 2 -PCR-A) for 84 genes related to the identification, growth, maintenance, and 
differentiation of SCs. 

Results: Cells from the corneal epithelium SC niche showed significant expression of ATP -binding cassette sub-family 
G member 2 (ABCG2) and cytokeratin (KRT)15, KRT14, and KRT5 genes. RT 2 -PCR-A results indicated an increased or 
decreased expression in 21 and 24 genes, respectively, in cells from the corneal SC niche compared to cells from the 
central corneal epithelium. Functional analysis by proprietary software found 4 different associated pathways and a novel 
network with the highest upregulated genes in the corneal SC niche. This led to the identification of specific molecules, 
chemokine (C-X-C motif) ligand 12 (CXCL12), islet- 1 transcription factor LIM/homeodomain (ISL1), collagen-type II 
alpha 1 (COL2A), neural cell adhesion molecule 1 (NCAM1), aggrecan (AC AN), forkheadbox A2 (FOXA2), Gap junction 
protein beta l/connexin32 (GJB1/Cnx32), and Msh homeobox 1 (MSX1), that could be used to recognize putative corneal 
epithelial SCs grown in culture and intended for transplantation. Other molecules, NCAM1 and GJB1/Cnx32, potentially 
could be used to positively purify them, and Par-6 partitioning defective 6 homolog alpha (PARD6A) to negatively purify 
them. 

Conclusions: Knowledge of these gene and molecular pathways has provided a better understanding of the signaling 
molecular pathways associated with progenitor-rich limbal epithelium. This knowledge potentially could give support to 
the design and development of innovative therapies with the potential to reverse corneal blindness arising from ocular 
surface failure. 



The cornea is the clear front of the eye through which 
light enters on its way to the retina. The corneal outer surface 
is covered by a stratified squamous nonkeratinized epithelium 
that resists constant attrition caused by exposure-induced 
dryness and potential light- induced damage [1]. To cope with 
this demand, constant renewal and maintenance of the corneal 
epithelium is achieved by stem cells (SCs) located at the 
circular border of the cornea in a region known as the 
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corneoscleral limbus. The basal epithelial cells of the limbal 
region are not homogeneous, but rather consist of diverse 
populations of SCs, transient amplifying cells, and terminally 
differentiated cells for which the total number and distribution 
are unknown [1-4]. Limbal SC deficiency (LSCD) syndrome 
occurs if limbal epithelial SCs (LESCs) are critically reduced 
and/or dysfunctional due to a multitude of conditions 
including genetic disorders (i.e., anirida), cicatrizing- 
autoimmune pathologies (i.e., Steven- Johnson syndrome, 
mucous membrane pemphygoid), severe infections, or 
external factors such as chemical or thermal burns, ultraviolet 
and ionizing radiation, contact lens wear, and multiple 
surgeries. The consequence of LSCD is a chronic pain 
inflammatory syndrome and loss of vision, greatly affecting 
quality of life and productivity [5]. 

Current treatment of LSCD relies on the inhibition of 
inflammation, protection, and provision of LESCs for 
reconstruction of the damaged corneas [5-7]. Strategies based 
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on transplantation of ex vivo expanded LESCs are becoming 
widely accepted today. The most frequently chosen technique 
includes harvesting autologous or allogenic limbal tissue that 
is then cultivated on amniotic membranes or fibrin matrices. 
Transplantation of these cultured cells has shown promising 
results [8-12]. However, it is usually not known what 
percentage of the transplanted cells is actually composed of 
SCs. It is likely that the success of each transplantation 
depends upon the number of SCs included. For example, 
enrichment of transplants with LESCs expressing the marker 
p63 increases the success rate [10]. It is therefore essential to 
improve the purity of the LESCs being transplanted to ensure 
good long-term transplantation results. 

Identifying LESCs is crucial for enrichment and 
characterization. Unfortunately, to date, no direct methods 
have been established because no single specific LESC 
marker is known. A variety of SC markers has been proposed 
to identify the LESC population. In addition, a diversity of 
differentiation markers has also been proposed to differentiate 
LESCs from terminally differentiated corneal epithelial cells 
[13-16]. Until now, the combination of positive and negative 
SC markers seems to be the most trustworthy way to 
characterize the putative SCs in the limbal epithelium. 
Typically, the major positive markers used are the 
transcription factor p63, the drug-resistance transporter ATP- 
binding cassette sub-family G member 2 (ABCG2), and some 
cytokeratins (KRTs) like KRT15 and KRT14. Among the 
most used as negative markers are KRT3 and KRT12, and the 
gap junction protein connexin 43, which are all typical of 
terminally differentiated cells [10,13,15,16]. 

Recently, great efforts have been made toward the 
identification of new molecular markers that may better 
distinguish LESCs from transient amplifying cells and 
terminally differentiated cells [16,17]. However, the variety 
of putative LESC markers and their role for the identification 
the LESC population is controversial [15,18]. The finding of 
new molecules that specifically identify LESCs would 
significantly enhance the purity of LESCs grown in culture 
and intended for transplantation. In addition, a better 
understanding of the molecular signaling pathways associated 
with the sternness of the limbal epithelium could facilitate a 
better diagnosis of LSCD and could also give support to the 
design and development of new and promising treatments. 
Therefore to discover new putative LESC markers, we 
analyzed the expression of 84 genes related to the 
identification, growth, maintenance, and differentiation of 
human SCs. Using a real time reverse transcription 
polymerase chain reaction array (RT 2 -PCR-A) with human 
corneal and limbal samples, we found increased and decreased 
expression of selected genes operating in 4 different pathways 
constituting signaling networks in the cells from the limbal 
stem cell niche. 
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METHODS 

Epithelial cell collection: Human tissue was used in 
accordance with the Declaration of Helsinki. Normal human 
corneoscleral tissues (n=24) were obtained 3 to 5 days post- 
mortem from the Barraquer Eye Bank (Barcelona, Spain). 
Limbal and central cornea epithelial cells were obtained using 
a modification of a previously described method [19-23]. In 
brief, a 7.5 mm trephine was used to isolate the cornea from 
the limbus, and the epithelium in the central button of the 
cornea was scraped to harvest differentiating epithelial cells 
for analysis of gene expression. Later, each corneoscleral rim 
was trimmed, and the endothelial layer and iris remnants were 
removed. The limbal rim was incubated with dispase II (5 mg/ 
ml; STEMCELL Technologies, Grenoble, France) at 37 °C 
for 2 h. The limbal epithelial sheets were then collected and 
treated with 0.25% trypsin with 0.03% EDTA (Invitrogen- 
Gibco, Inchinnan, UK) at 37 °C for 10 min to isolate single 
cells. There were, therefore, 24 samples of 2 different types 
of epithelial cells: differentiating corneal epithelial cells and 
stem cell-containing population of limbal epithelial cells 
derived from the corneal epithelial stem cell niche. 
RNA isolation and reverse transcription: Total RNA was 
extracted by Qiagen RNeasy Mini Kit (QIAGEN Inc., 
Valencia, CA) under standard conditions, and treated with 
RNase-free DNase following our previously described 
method [24-26]. Briefly, samples were collected in RNA lysis 
buffer (1:100 (3-mercaptoethanol-buffer RLT), purified in 
QIAshredder columns, and treated with RNase-Free DNase I 
Set (QIAGEN Inc.) following the manufacturer's instructions. 
Agarose gel electrophoresis and ethidium bromide staining 
were used to check the integrity and size distribution of the 
purified RNA. The first strand of cDNA was synthesized with 
random hexamer using M-MuLV Reverse Transcriptase 
(Amersham Pharmacia Biotech Europe GmbH, Barcelona, 
Spain) [24-26]. 

Real time polymerase chain reaction (rt-PCR): The cDNA 
from the limbal and corneal epithelial cells was mixed with 
Taqman assay primers and minor groove binder probes 
specific for glyceraldehyde 3 -phosphate dehydrogenase 
(GAPDH), KRT3, KRT5, KRT7, KRT12, KRT14, KRT15, 
KRT19, p63 and ABCG2 (Table 1) and with a Taqman 
Universal PCR Master Mix AmpErase UNG (Applied 
Biosystems, Foster City, CA) in a 7500 Real Time PCR 
System (Applied Biosystems) according to the previously 
described method [27-31]. An aliquot of 2 ul containing 20 
ng of cDNA was used for PCR in a total volume of 20 ul 
containing: 7 ul double-distilled water, 1 ul of 20x target 
primers and probe, 1 0 ul of 2* Taqman Universal PCR Master 
Mix. PCR parameters consisted of uracil N-glycosylane 
activation at 50 °C for 2 min, pre-denaturation at 95 °C for 2 
min, followed by 40 cycles of denaturation at 95 °C for 15 s, 
and annealing and extension at 60 °C for 1 min. 

Assays were performed in triplicate. A nontemplate 
control and total RNA without retrotranscription were 
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Table 1. Oligonucleotide primers and probes used for real time PCR 
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^Identification number from Applied Biosystems (www.appliedbiosystems.com). 



included in all experiments to evaluate PCR and DNA 
contamination of the reagents. GAPDH was used as an 
endogenous reference for each reaction to correct for 
differences in the amount of total RNA added. To verify the 
validity of using GAPDH as an internal standard control, the 
efficiencies of the genes and GAPDH amplifications were 
compared. 

The comparative cycle threshold (Ct) method, where the 
target fold=2~ AACt , was used for analyzing the results (Applied 
Biosystems User Bulletin, No.2, P/N 4303859) [27-31]. 
Corneal mRNA served as the calibrator control. The results 
were reported as a fold upregulation when the fold-change for 
limbal cells was greater than one compared to corneal cells. 
If the fold-change was less than one, the negative inverse of 
the result was reported as a fold down-regulation. Significant 
differences (p<0.05) were evaluated by Student's Mest. 
Real time PCR array: The samples were pooled, creating 4 
groups of 6 each, and used for further study. Analysis using a 
real time PCR (rt-PCR) array was performed according to the 
manufacturer's recommendations using the Human Stem Cell 
RT 2 Profiler™ (Super Array Bioscience, Izasa, S.A., 
Barcelona, Spain) that used SYBR® Green I dye detection. 

We studied the expression of 5 housekeeping genes, 3 
RNAs and PCR quality controls, and 84 human genes related 
to: 

1. SC specific markers (cell cycle regulators, 
chromosome and chromatin modulators, genes 
regulating symmetric/asymmetric cell division, self- 
renewal, cytokines and growth factors, genes 
regulating cell-cell communication, cell adhesion 
molecules and metabolism), 

2. SC differentiation markers (embryonic, 
hematopoietic, mesenchymal, and neural cell lineage 
markers), and 

3. Signaling pathways important for SC 
maintenance (Notch and Wnt pathways). 



The following components were mixed in a 5 -ml tube: 
1,275 |il of the 2x Super Array PCR Master Mix, 102 ul (100 
ng) of the diluted first strand cDNA synthesis reaction, and 
1,173 ul double-distilled H2O. This mixture and template 
cocktail (25 ul each) was added to each well of the PCR array. 
Real time PCR (7500 Real Time PCR System) was then 
performed as follows: 10 min at 95 °C, 40 cycles of 15 s at 
95 °C, and 1 min at 60 °C. Assays were performed in duplicate. 
A melting curve program was run and a dissociation curve 
was generated for each well in the entire plate to verify the 
identity of each gene amplification product. 

For data analysis, the Ct method was performed using an 
Excel-based PCR Array Data Analysis template that was 
downloaded from the SuperArray website. This program 
automatically performed the following calculations and 
interpretation of gene expression based upon threshold cycle 
data from a real-time instrument: 

1. Changed to 35 all Ct values greater than 35 and 
Ct values not detected. At this point, any Ct value 
equaled to 35 was considered a negative call. 

2. Examined the threshold cycle values of the 
genomic DNA control, reverse transcription control, 
and positive PCR control wells. 

3. Calculated the ACt for each gene in each plate. 

We used the average of the five housekeeping gene Ct 
values as a normalization factor. The results are reported as a 
fold upregulation or down-regulation in the same way as 
previously explained for real time PCR (above). 

Pathway analysis: Excel spreadsheets containing gene 
identifier lists together with the corresponding expression 
values were uploaded into Ingenuity Pathways Analysis 
(IPA; Ingenuity® Systems, Redwood City, CA) to identify 
relationships among the genes of interest. The basis of the IPA 
program consisted of the Ingenuity Pathways Knowledge 
Base (IPKB) that was derived from known functions and 
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Epithelial cells from the limbal stem cell niche vs central corneal cells 
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Figure 1. Relative expression of 
cytokeratins (KRTs), p63, and ABCG2 
genes. Real time PCR was performed in 
24 samples of limbal and corneal 
epithelial cells. Mean mRNA 
expression in limbal epithelial cells was 
expressed relative to that corneal 
epithelial cells. Positive values 
indicated relatively greater expression 
in limbal cells, while negative values 
indicated relatively less expression in 
limbal. Significance was analyzed by 
Student's f-test. *p<0.05, **p<0.001, 
***p<0.00001. 



interactions of genes published in the literature. Thus, the IPA 
tool allowed the identification of biologic networks, global 
functions, and functional pathway(s) of a particular data set. 
Each gene identifier was mapped to its corresponding gene 
object in the IPKB. Networks of the genes were then 
algorithmically generated based on their connectivity. 

Each gene product was assigned to functional and sub- 
functional categories. IPA software then used the associated 
library of canonical pathways to identify the most significant 
ones in the data set. Benjamini-Hochberg multiple testing 
correction was used to calculate a p-value to determine the 
probability that each biologic function or canonical pathway 
assigned to the data set was due to chance alone. In addition, 
significance of the association between the data set and the 
canonical pathway was calculated as a ratio of the number of 
genes from the data set that mapped to the pathway divided 
by the total number of genes that map to the canonical 
pathway. The 'Pathway Designer' tool of the IPA software 
was used for the graphical representation of the molecular 
relationships between gene products. Gene products were 
represented as nodes, and the biologic relationship between 
two nodes was represented as an edge (line). All edges were 
supported by at least one reference from the literature, from a 
textbook, or from canonical information stored in the IPKB. 

RESULTS 

Real time PCR analysis for corneal and limbal epithelial cell 
markers: To select the purest population of corneal and limbal 
epithelial cells, we performed rt-PCR assays to evaluate the 
expression of markers considered to be abundant in the limbal 
stem cell niche. These markers included KRT14, KRT15, 
ABCG2, and transcription factor p63 [13,20,32-35]. For 
terminally differentiated corneal epithelial cells, we looked 
for the expression of KRT3,KRT7, and KRT1 2 [3,36], as well 
as for other cytokeratins like KRTS and KRT19 [ 1 5] . In the 24 
samples analyzed, all of the studied KRT genes were 
expressed (Figure 1). In the limbal epithelial cells, expression 
was significantly reduced for most cytokeratin genes that are 



normally expressed in large amounts in terminally 
differentiated epithelial cells [15,35]. The reductions for 
KRT3, KRT 7, KRT12, mdKRT19, which varied between 2.03 
and 3.54 fold, were all significant except for KRT 12 (p<0.05 
for KRT3 and KRT1 7, p>0.05 for KRT1 2, andpO.00001 for 
KRT19, Figure 1). In contrast, KRT5, KR Tl 4, and KRT15 were 
more highly expressed in the limbal than the corneal epithelial 
samples, with increases ranging from 2.29 to 29.46 fold 
(p<0.05, O.001, and O.00001, respectively, Figure 1). 

Gene expression of associated LESC niche markers 
ABCG2 and p63 were found in all of the samples analyzed. 
Expression levels of ABCG2 were 39.1 fold greater in the 
limbal epithelial cell population than in the corneal epithelial 
one (p<0.00001, Figure 1). However, expression of 
transcription factor p63 was the same in both cell populations. 

Summarizing our results so far, the purest SC-containing 
population of limbal epithelial cells had significantly higher 
expression of ABCG2 (39 fold), KRT 15 (29.5 fold), KRT 14 
(5.6 fold), and KRT5 (2.3 fold) than did the corneal epithelial 
cell population. Furthermore, the limbal cells had 
significantly lower expression of KRT3, KRT7, and KRT19. 
Neither KRT 12 nor p63 were useful as gene markers to 
differentiate between the two cell populations. 

Real time PCR array: The 24 samples previously analyzed by 
real time PCR were pooled to perform the PCR array. The 
dissociation curve was analyzed for the 84 genes studied, and 
no DNA contamination was detected. The results indicated 
increased expression of 21 genes and decreased expression of 
24 genes for limbal cells compared to corneal cells. Eleven 
genes had a greater than ninefold increased expression and 10 
genes had a greater than fourfold decreased expression (Table 
2). 

Among the 1 1 most upregulated genes (Table 2) from the 
limbal SC niche, three coded for extracellular space proteins 
(chemokine [C-X-C motif] ligand 12 [CXCL12], collagen- 
type II alpha 1 [COL2A], and aggrecan [ACAN]), three for 
transcription factors located at the nucleus (islet- 1 
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Table 2. Up- and down-regulated genes in cells of the lumbal epithelial stem cell niche. 
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Genes with higher (+) and lesser (-) expression in limbal epithelial cells compared to terminally differentiated corneal epithelial 
cells. Fold change was calculated by PCR array using the comparative Ct method. ^Indicates cellular location where protein is 
expressed. 



transcription factor LIM/homeodomain [ISL1], forkhead box 
A2 [FOXA2], and Msh homeobox 1 [MSX1]), four for plasma 
membrane proteins (neural cell adhesion molecule 1 
[NCAM1], ABCG2, Gap junction protein beta 1 [GJB1], and 
CD 8b molecule [CD8B]), and only one for a cytoplasmic 
protein (KRT15). Among them, the most upregulated 
expression was for the chemokine CXCL12 gene with 26.45 
fold increased expression. 

Among the 10 most down-regulated genes (Table 2), two 
coded for extracellular space proteins (Desert hedgehog 
homolog [DHH] and Jagged 1 [JAG1]), two for nuclear 
proteins (Cell division cycle 2 [CDC2] and Cyclin A2 
[CDCNA2]), two for plasma membrane proteins (Gap 
junction protein alpha 1 [GJA1] and Par-6 partitioning 
defective 6 homolog alpha [PARD6A]), and four for 
cytoplasmic proteins (K[lysine] acetyltransferase 2A 
[KAT2A], Dishevelled dsh homolog 1 [DVL1], S100 calcium 
binding protein B [SIOOB], and Frequently rearranged in 



advanced T-cell lymphomas [FRAT1]). Among them, the 
most down-regulated expression was for the DHH peptidase 
gene with 17.66 fold decreased expression. 

Signaling pathways — Seventy canonical signaling 
pathways were significantly affected across the entire data set 
identified by IP A (Table 3, Figure 2). The highest upregulated 
gene was SOX (9.2 fold, Figure 2) in the Wnt/p-catenin 
signaling pathway, also known as SRY (sex determining 
region Y)-box 2, Entrez Gene 6736). The most down- 
regulated gene was GJA1 (6.9 fold), also known as gap 
junction protein, alpha 1 (Entrez Gene 2697). 

Predicted functional effects — The IPA program 
determined if groups of genes with significantly changed 
expression levels were associated with altered biologic 
functions and diseases (Table 4). Here IPA identified 71 
functional categories that were significantly affected. The 
most prominent cellular and molecular functions implicated 
were cellular development, cell death, gene expression, 
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Table 3. Ingenuity canonical pathways that were most significantly affected. 



Ingenuity Canonical Pathways 
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Sonic Hedgehog Signaling 
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Calcium Signaling 
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CCR5 Signaling in Macrophages 

Neurotrophin/TRK Signaling 

Caveolar-mediated Endocytosis 
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LPS-stimulated MAPK Signaling 

Bile Acid Biosynthesis 

Chemokine Signaling 



-Log(B-H P-value)* 

8.54E+00 



3.28E+00 

2.58E+00 

1.49E+00 

1.34E+00 

1.34E+00 

1.34E+00 

1.27E+00 
1.27E+00 

1.11E+00 
1.07E+00 

8.86E-01 
8.86E-01 
8.21E-01 
7.08E-01 
7.08E-01 
7.08E-01 
6.12E-01 
5.94E-01 

4.92E-01 
4.92E-01 

4.92E-01 
4.73E-01 
4.34E-01 
4.12E-01 
3.39E-01 
3.39E-01 
3.39E-01 
3.39E-01 
3.39E-01 
3.39E-01 
3.39E-01 

3.39E-01 
3.39E-01 
3.39E-01 
3.39E-01 
3.39E-01 
3.39E-01 
3.39E-01 
3.39E-01 



Ratio 

6.67E-02 



9.76E-02 

6.78E-02 

2.20E-02 

2.58E-02 

2.42E-02 

1.52E-02 

3.49E-02 
2.07E-02 

2.78E-02 
4.65E-02 

L83E-02 
2.08E-02 
3.28E-02 
1.55E-02 
2.50E-02 
2.17E-02 
2.02E-02 
1.92E-02 

1.65E-02 
3.85E-02 

1.09E-02 
1.48E-02 
3.23E-02 
1.22E-02 
9.71E-03 
1.56E-02 
1.01E-02 
8.62E-03 
9.01E-03 
1.54E-02 
L37E-02 

1.43E-02 
1.16E-02 
1.32E-02 
1.25E-02 
1.16E-02 
1.27E-02 
1.03E-02 
1.30E-02 



Molecules 

SOX2, CDH2, GJA1, 
AXIN1,FRAT1,DVL1, 
BTRC, FZD1,CCND1, 

WNT1, APC 
NOTCH2, DLL1, DTX1, 
JAG1 

CCNE1,HDAC2, BTRC, 

CCND1 
FGF4, CDC42, ACTC1, 
APC, FGF1 
CCNA2, CCNE1, 
ALDH1A1, CCND1 
FGF4, CDC42, ACTC1, 
FGF1 

CDC42, BMP2, FZD1, 

WNT1,BMP1 
FGF4, FGFR1, FGF1 
CDC42, AXIN1, CXCL12, 
FGF1 
CD8A, CD3D, CD8B 
BTRC, CDC2 

CDC42, ACTC1, PARD6A 
HDAC2, BMP2, BTRC 
HDAC2, CD3D 
CDC42, CXCL12, ACTC1 
BMP2, BMP1 
CDC42, ACTC1 
CDC42, CCND1 
CDC42, ACTC1 

CDC42, CD3D 
CD3D 

HSPA9, BGLAP, CD3D 
FGFR1, FGF1 
CDC2 
ALDH1A1 
HDAC2, ACTC1 

ALDH1A1 
CDC42, ACTC1 
HDAC2, HSPA9 
ALDH1A1 
CCND1 
ADAR 

CDC42 

CD3D 

CDC42 

ACTC1 
ALDH1A1 

CDC42 
ALDH1A1 
CXCL12 
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Table 3. Continued. 



Ingenuity Canonical Pathways 

L 2 -alanine Metabolism 
VDR/RXR Activation 
Butanoate Metabolism 
Pyruvate Metabolism 

CTLA4 Signaling in Cytotoxic T Lymphocytes 

TGF-P 2 Signaling 

Lysine Degradation 

Propanoate Metabolism 

Apoptosis Signaling 

p53 Signaling 

VEGF Signaling 

SAPK/JNK Signaling 

Glycerolipid Metabolism 

Valine, Leucine and Isoleucine Degradation 

FXR/RXR Activation 

Glycolysis/Gluconeogenesis 

Nicotinate and Nicotinamide Metabolism 

fMLP Signaling in Neutrophils 

Arginine and Proline Metabolism 

PI3K/AKT Signaling 

Protein Ubiquitination Pathway 

Inositol Phosphate Metabolism 

B Cell Receptor Signaling 

CXCR4 Signaling 

IL-8 Signaling 

Tryptophan Metabolism 

RAR Activation 

Role of NFAT in Regulation of the Immune 
Response 

Fatty Acid Metabolism 

NRF2-mediated Oxidative Stress Response 

LPS/IL-1 Mediated Inhibition of RXR Function 



-Log(B-H P 

3.39E- 
3.39E- 
3.39E- 
3.39E- 
3.39E- 
3.39E- 
3.39E- 
3.39E- 
3.39E- 
3.39E- 
3.39E- 
3.39E- 
3.39E- 
3.39E- 
3.39E- 
3.33E- 
3.33E- 
3.03E- 
3.01E 
3.00E- 
2.60E- 
2.60E- 
2.54E- 
2.51E- 
2.51E- 
2.51E- 
2.51E- 
2.51E- 



-value)* 

01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 



2.51E-01 
2.51E-01 
2.29E-01 



Ratio 

1.01E-02 
1.25E-02 
7.75E-03 
6.90E-03 
1.16E-02 
1.16E-02 
6.94E-03 
7.94E-03 
1.06E-02 
1.15E-02 
1.05E-02 
1.08E-02 
6.90E-03 
9.35E-03 
L00E-02 
7.09E-03 
7.75E-03 
8.00E-03 
5.62E-03 
7.41E-03 
4.95E-03 
5.78E-03 
6.54E-03 
6.10E-03 
5.46E-03 
4.20E-03 
5.49E-03 
5.38E-03 

5.29E-03 
5.46E-03 
5.05E-03 



Molecules 

ALDH1A1 

BGLAP 
ALDH1A1 
ALDH1A1 
CD3D 
BMP2 
ALDH1A1 
ALDH1A1 
CDC2 
CCND1 
ACTC1 
CDC42 
ALDH1A1 
ALDH1A1 

FOXA2 
ALDH1A1 
CDC2 
CDC42 
ALDH1A1 
CCND1 
BTRC 
CDC2 
CDC42 
CXCL12 
CCND1 
ALDH1A1 
ALDH1A1 
CD3D 

ALDH1A1 

ACTC1 
ALDH1A1 



*The p-value was calculated using the Benjamini-Hochberg (BH) method. The ratio was calculated as the number of genes in 
a given pathway that met the cutoff criteria, two in this case, divided by the total number of genes that made up the pathway. 



cellular assembly and organization, and cellular growth and 
proliferation. The most frequent significant physiologic 
system developments were tissue, organismal, embryonic, 
nervous system, and organ development. 

Gene networks — The IPA program constructed 4 gene 
networks that were significantly interconnected. The first 
network (Figure 3A) contained 22 genes concerned with the 
auditory and vestibular system development and function, 
organ development, and cancer. Upregulated genes included 
FGF4, FGFR1, ISL1, MSX1, NCAM1, NOTCH2, SOX2, T, 
and WNTL Down-regulated genes included APC, AXIN1, 
BGLAP, BTRC, CCND1, CDH2, DLL1, DTX1, DVL1, 
FGF1, FRAT1, HSPA9, mdJAGL 

The second network (Figure 3B) contained 17 genes 
associated with cancer, connective tissue development and 
function, and skeletal and muscular system development and 
function. Upregulated genes included ACAN, BMP2, CD8B, 



COL2A1, and CXCL12. Down-regulated genes included 
ACTC1, BMP1, CCNA2, CCNE1, CD3D, CD8A, CDC42, 
FZD1, GJA1, GJB2, PARD6A, and SIOOB. 

ALDH1 occupied a central position in the third network 
(Figure 4A), which contained 12 genes concerned with drug 
metabolism, small molecule biochemistry, and cell 
morphology expression. Upregulated genes included 
ABCG2, COL9A1, FOXA2, KRT15, and NEUROG2. Down- 
regulated genes included ADAR, ALDH1A1, ASCL2, 
COL9A2DHH, GDF3, GJB2, and OPRSL 

Finally, the fourth network (Figure 4B) contained 6 genes 
affecting cancer, cell cycle, and skeletal and muscular 
disorders. There were 2 upregulated genes, GJB1 and MME, 
and 4 down-regulated, CDC2, HDAC2, KAT2A, and MYST2. 

Customized gene network — Using the IPKB, we 
explored possible functional relationships among the six 
highest upregulated limbal epithelium progenitor-rich cell 
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6.8 




Figure 2. Wnt/p-catenin signaling pathway generated by Ingenuity Pathway Analysis (IP A). The IP A depicted the genes involved, their 
interactions, and the cellular and metabolic reactions that constituted the pathway. Colored molecules represented genes that appeared in the 
data set studied. Red and green molecules were up- and down-regulated, respectively, in limbal epithelial cells. Gray molecules did not meet 
the user defined cutoff of 2. 



genes (Table 2): (1) CXCL12, (2) ISL1, (3) COL2A1, (4) 
NCAM1, (5) ABCG2, and (6) KRT15. We obtained a network 
with 29 genes. The protein products of 14 genes were active 
in the nucleus, one in the cytoplasm, six in the plasma 
membrane, and seven in the extracellular space (Figure 5). 

CXCL12, also called stromal cell-derived factor 1 
(SDF1), encodes for small cytokines that belong to the 
intercrine family (Entrez Gene 6387). We chose it as the 
central gene in the network because in humans it directly or 
indirectly interacts with the other genes that we added. Among 
the 6 most upregulated genes, only ISL1, which encodes for a 
member of the LIM/homeodomain family of transcription 



factors and may play an important role in regulating insulin 
gene expression (Entrez Gene 3670), did not have any 
connections with other genes in this network. 

DISCUSSION 

Isolation and characterization of tissue specific SCs to study 
their functional properties is one of the main research 
aspirations for regenerative medicine. In the context of ocular 
surface therapy, the ability to identify, purify, and characterize 
LESCs is an essential goal. However, the lack of LESC 
specific markers has been an obstacle for their isolation and 
subsequent biologic and functional characterization. Using 
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Table 4. Molecules significantly associated with relevant functions and diseases. 



Category B-H P-value Molecules 

Cellular Development 3.49E-13-1.53E-02 37 

Tissue Development 3.61E-10-1.53E-02 31 

Organismal Development 3.91E-09-4.75E-03 22 

Embryonic Development 2.37E-08-1.53E-02 26 

Cell Death 3.61E-08-1.53E-02 35 

Gene Expression 2.67E-07-1.53E-02 32 

Cellular Assembly and Organization 2.67E-07-1 .53E-02 22 

Nervous System Development and Function 2.67E-07-1.53E-02 29 

Cancer 3.1E-07-1.53E-02 40 

Cellular Growth and Proliferation 5.45E-07-1 .53E-02 39 

Cell Morphology 5.45E-07-1.53E-02 23 

Cell Cycle 1.15E-06-1.53E-02 17 

Organ Development 1 .69E-06-1 .36E-02 20 

Skeletal and Muscular Disorders 3.01E-06-1 .53E-02 1 1 

Renal and Urological Disease 4.17E-06-1.36E-02 7 

Genetic Disorder 5.66E-06-1.53E-02 39 

Developmental Disorder 7.55E-06-1 .53E-02 1 7 

Connective Tissue Development and Function 4.09E-05-1.53E-02 21 

Skeletal and Muscular System Development and Function 4.09E-05-1 .5E-02 1 7 

Hematological System Development and Function 8.41E-05-1 .53E-02 1 8 

Hematopoiesis 8.41E-05-1.53E-02 13 

Neurologic Disease 1 .35E-04-1 .53E-02 27 

Lymphoid Tissue Structure and Development 1 .86E-04-1 .53E-02 9 

Auditory and Vestibular System Development and 2.5 1E-04-8.77E-03 5 
Function 

Cell-To-Cell Signaling and Interaction 2.92E-04-1 .53E-02 23 

Cellular Movement " 2.98E-04-1.53E-02 22 

Cellular Function and Maintenance 3.24E-04-1 .53E-02 14 

Cardiovascular System Development and Function 3.82E-04-1 .53E-02 14 

Hepatic System Disease 3.82E-04-1.53E-02 14 

Reproductive System Disease 5.45E-04-1.53E-02 19 

Gastrointestinal Disease 5.48E-04-1.53E-02 21 

Tissue Morphology 6.35E-04-1 .53E-02 1 8 

Energy Production 6.56E-04-6.56E-04 3 

Molecular Transport 6.56E-04-1.53E-02 3 

Nucleic Acid Metabolism 6.56E-04-1 .53E-02 4 

Small Molecule Biochemistry 6.56E-04-1 .53E-02 9 

Organ Morphology 1 .04E-03-1 .53E-02 1 3 

Tumor Morphology 1.13E-03-1.53E-02 10 

Metabolic Disease 1 .33E-03-1 .33E-03 3 

DNA Replication. Recombination, and Repair 1 .5E-03-1 .53E-02 1 1 

Connective Tissue Disorders 1.54E-03-7.79E-03 10 

Humoral Immune Response 2.21E-03-3.6E-03 4 

Visual System Development and Function 2.21E-03-7.79E-03 4 

Psychological Disorders 2.25E-03-2.25E-03 6 

Infection Mechanism 2.47E-03-9.23E-03 4 

Post-Translational Modification 2.69E-03-1 .53E-02 1 1 

Carbohydrate Metabolism 2.81E-03-2.81E-03 5 

Lipid Metabolism 2.81E-03-2.81E-03 3 
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Table 4. Continued. 



Category 


B-H P-value 


Molecules 


Drug Metabolism 


2.81E-03-1.53E-02 


3 


Endocrine System Development and Function 


2.81E-03-1.53E-02 


4 


Hair and Skin Development and Function 


3.04E-03-1.53E-02 


6 


Reproductive System Development and Function 


3.45E-03-1.53E-02 


4 


Hematological Disease 


5.19E-03-1.53E-02 


8 


Cardiovascular Disease 


5.28E-03-5.28E-03 


4 


Cell-mediated Immune Response 


6.31E-03-1.3E-02 


5 


Organismal Injury and Abnormalities 


6.88E-03-6.88E-03 


3 


Digestive System Development and Function 


7.14E-03-7.14E-03 


3 


Organismal Survival 


7.53E-03-7.53E-03 


12 


Hepatic System Development and Function 


8.76E-03-8.76E-03 


2 


Respiratory System Development and Function 


8.76E-03-1.53E-02 


3 


Antigen Presentation 


1.23E-02-1.23E-02 


2 


Inflammatory Disease 


1.53E-02-1.53E-02 


2 


Cell Signaling 


1.53E-02-1.53E-02 


2 


Protein Trafficking 


1.53E-02-1.53E-02 


2 


Vitamin and Mineral Metabolism 


1.53E-02-1.53E-02 


2 


Renal and Urological System Development and Function 


1.53E-02-1.53E-02 


2 


Auditory Disease 


1.53E-02-1.53E-02 




Cellular Compromise 


1.53E-02-1.53E-02 




Dermatological Diseases and Conditions 


1.53E-02-1.53E-02 




Immune Cell Trafficking 


1.53E-02-1.53E-02 




RNA Post- Transcriptional Modification 


1.53E-02-1.53E-02 





*The p-value was calculated using the Benjamini-Hochberg (BH) method. 



cells isolated from the limbal SC niche, we compared the 
expression profile of 84 SC phenotype-related genes with 
cells from the differentiating central corneal epithelium zone. 
Our goal was to provide new information for molecules that 
are predominantly expressed in the stem cell-containing 
population of human limbal epithelial cells. Knowledge 
regarding these LESC potential markers could be used to 
enhance isolation of the cells and develop a better 
understanding of their biologic functions. 

To know the gene expression pattern of the isolated cell 
samples, we first performed a PCR analysis for corneal and 
limbal epithelial markers. The limbal epithelial cell 
population expressed high levels of ABCG2, KRT5, KRT14 
and KRT15 and low levels of KRT3, KRT7, and KRT19. 
Unexpectedly, we did not find significant differences between 
limbal and corneal epithelial cells for transcription factor 
p63 expression. In 2001, Pellegrini et al. [33] proposed p63 
as the first positive marker of LESCs. This has generated a 
certain level of controversy because several groups have since 
found that p63 is also expressed by most of the terminally 
differentiated basal epithelial cells throughout the cornea 
[18,37,38]. Our findings are consistent with the idea that p63 
is not specific enough to be a definitive marker for LESCs, 
although perhaps it could be helpful for identifying 



incompletely differentiated corneal epithelial cells [18]. It is 
worth noting that the a isoform of ANp63 has been proposed 
to be a rather more specific and useful marker for LESCs than 
the other isoforms of this transcription factor [10,39]. 

Several microarray studies have attempted to identify 
markers and signaling pathways associated with different 
ocular surface cell phenotypes [32,40-49]. We chose the 
RT 2 -PCR-A system because it utilizes real-time PCR in 
combination with microarray analysis to detect the 
simultaneous expression of many genes. We used IPA to 
analyze our results from the PCR array, creating three 
different analysis types that responded to three different 
questions: (1) What well characterized cell signaling and 
metabolic canonical pathways are most relevant to our data 
set? (2) What regulatory networks exist among the genes and 
proteins of our data set? (3) What previously unknown, unique 
customized networks that can serve as biologic models are 
present in our data set? 

Among the 84 genes we studied, 11 were highly 
upregulated and 10 were highly down-regulated; however less 
highly regulated genes may also be important in relation to 
SC properties. The most highly expressed in the limbal 
epithelium progenitor-rich cells compared to central corneal 
epithelial cells was the chemokine CXCL2. To explore 
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Figure 3. Networks generated by IPA related to the development and the function of the auditory, vestibular, skeletal and muscular systems 
and to the cancer development. Auditory and vestibular system development and function, organ development, and cancer network (A), and 
cancer, connective tissue development and function, skeletal and muscular system development and function network (B) generated by IPA. 
The networks contained nodes composed of genes/gene products and edges that indicated a relationship between the nodes in the cellular and 
subcellular locations indicated. Classes of nodes were indicated by shape to represent different functionalities. Colored molecules represented 
genes that appeared in the data set studied. Red and green molecules were upregulated and down-regulated, respectively, in the limbal epithelial 
cells. Gray molecules did not meet the user defined cutoff of 2. White indicated the molecule was added from the IPKB. 



molecular signatures of progenitor cells, we further analyzed 
six highly expressed genes, starting with the chemokine 
CXCL12, to create our customized gene network with a total 
of 29 molecules. 

Chemokines are 8- to 10-kDa proteins that are potent 
activators and chemoattractants for different leukocyte 
subpopulations and some non-hematopoietic cells such as 
epithelial cells, fibroblasts, and endothelial cells [50]. 
CXCL12 and its receptor CXCR4 are expressed in cultured 
human corneal fibroblasts [51]. They may play a key role in 
angiogenesis and be involved in ocular neovascularization as 
well as in the recruitment of inflammatory or vascular 
endothelial cells to sites of corneal injury. In a recent 
microarray analysis of pig limbal side population cells, 
CXCR4 had the greatest overexpression ratio [42]. CXCR4 is 
also upregulated in pig and human conjunctiva side 
population cells [41,42]. Based on all of these findings, the 
CXCL12/CXCR4 pair could serve as a suitable marker to 
identify ocular surface SCs in a species-independent way. 
CXCL12/CXCR4 signaling is also critical for the 
mobilization and recruitment of mesenchymal SCs (MSCs) to 
infarcted hearts and fracture sites in bones [52,53]. 
Additionally, Ye et al. [54] recently reported that systemically 
transplanted bone marrow MSCs can engraft to injured cornea 
and promote wound healing by differentiation, proliferation, 
and synergizing with hematopoietic SCs. Thus we 



hypothesize that corneal homing of MSCs after ocular surface 
wounding could be mediated by release of CXCL12 from 
limbal epithelial cells and corneal fibroblast. Potentially, 
CXCL12 topical administration could be used to enhance 
MSC homing to injured corneal and limbal areas, facilitating 
the regenerative processes. 

In addition to locating the SCs of the epithelium, the ideal 
SC marker should also allow for isolation and enrichment of 
viable SCs from a heterogeneous epithelial cell population. 
For that reason, cell surface proteins such as cell-cell and cell- 
matrix adhesion molecules, as well as cell surface receptors, 
may be the best candidates for new positive and negative 
putative LESC markers. Based on our results and others [15, 
20], the plasma membrane transporter ABCG2 appears to be 
the most useful cell surface marker for the identification and 
isolation of LESCs. 

An example of a negative potential marker, one that 
indicates the absence of SC properties, could be PARD6A. 
This gene is a member of the PAR6 family and encodes a cell 
membrane protein involved in the control of epithelial cell 
polarity and tight junction assembly [55,56] and in epithelial- 
to-mesenchymal transition [57]. Expression of PARD6A in 
cells from the limbal stem cell niche was reduced fivefold 
compared to the corneal epithelial cells. 

Another such negative marker is the gap junction protein 
connexin 43 (GJA1) that is abundantly expressed in the 
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Figure 4. Networks generated by IPA related to drug metabolism, small molecule biochemistry and cell morphology and to cancer, cell cycle, 
skeletal and muscular disorders. Drug metabolism, small molecule biochemistry and cell morphology network (A), and cancer, cell cycle, 
skeletal and muscular disorders network (B) generated by IPA. The network contained nodes (gene/gene product) and edges (indicating a 
relationship between the nodes) showing the cellular/subcellular location as indicated. Function classes of nodes were indicated by shape to 
represent functional class. Colored molecules represented genes that appeared in the data set studied. Red and green molecules were upregulated 
and down-regulated, respectively, in limbal epithelial cells. Gray molecules did not meet the user defined cutoff of 2. White indicated the 
molecule was added from the IPKB. 



corneal but not in the limbal epithelium [19,58,59]. Membrane 
channel connexins (Cxs) form gap junctions that have been 
implicated in the homeostatic regulation of multicellular 
systems [60]. It is assumed that SCs of the limbal epithelium 
lack connexins and metabolite transfer capacity due to 
apparent self-sufficiency and absence of necessity for direct 
cell-to-cell communication [58]. However, our results showed 
upregulated expression of a related gene, Cx32 (GJB1), in 
limbal cells which was reported to be absent in human corneal 
epithelial cells [46]. Furthermore, Figueira et al. [32] recently 
described the expression of Cx32 in human fetal limbus and 
in cultured adult primary limbal explant epithelium. Similarly, 
hematopoietic cells were assumed not to express Cxs; 
however, hematopoietic SCs express Cx32 in response to 
chemical insult and also while maintaining the quiescent, 
noncycling state of primitive hematopoietic progenitor cells 
[61,62]. Although further investigations are required to 
confirm the role of Cx32 in LESCs, we propose this cellular 
surface protein as a new putative positive marker for the 
identification and isolation of human LESCs. 

Expression of the neural cell adhesion molecule 1 
(NCAM1) was highly upregulated in the limbal epithelial 
cells. NCAM is broadly expressed during development and 
plays a essential role in cell division, migration, and 
differentiation [63]. A decrease in NCAM expression during 



the development of the ocular lens has been associated with 
lens epithelial cell differentiation [64]. However NCAM is 
also expressed in cells of many fully developed tissues and 
organs including the cornea and lens epithelium [65]. For that 
reason, we believe it is not specific enough to serve as a 
potential single LESC marker. 

The limbal epithelium may contain a higher proportion 
of immune-related cells such as macrophages, lymphocytes, 
and antigen presenting cells than does the central corneal 
epithelium [66,67]. Thus the presence of significant portions 
of marker transcripts derived from these kinds of cells is not 
surprising. The best example of this is CD 8, a plasma 
membrane specific marker of T cells [68], that was 
overexpressed in the limbal-derived cells. This confirms the 
greater presence of immune-related cells in the limbal 
epithelium than in the corneal epithelium [66,67]. 

Analysis of our RT 2 -PCR-A data with IPA software 
recognized that the most significantly affected canonical 
pathway was Wnt/(32-catenin signaling, consistent with the 
recent findings of Bian et al. [43]. Wnt signaling is involved 
in practically every aspect of embryonic development and also 
controls homeostatic self-renewal in several adult tissues 
[69]. Among the studied molecules that belong to this 
pathway, SOX2 and Wnt were the highest upregulated genes, 
9.2 and 6.8 fold, respectively. The SOX2 gene encodes a 
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Figure 5. Customized gene network based upon the six most highly upregulated limbal epithelial cell genes. We explored possible functional 
relationships between the six highest upregulated limbal epithelial cells genes (in red) using the IPKB. Our customized pathway contained 
nodes composed of genes/gene products and edges that indicated a relationship between the nodes in the cellular and subcellular locations 
indicated. White indicates that the molecule was added from the IPKB. 



member of the SRY-related HMG-box (SOX) family of 
transcription factors implicated in the regulation of embryonic 
development and in the determination of cell fate [70]. Wnt 
signaling is required for the establishment of hair follicles, 
playing a key role in the activation of bulge SCs to progress 
toward hair formation [69,71]. Zhou et al. [48] prepared a 
transcriptional profile of mouse limbal and corneal epithelial 
basal cells. Consistent with our results, they found elevated 
expression of certain genes that were also upregulated in the 
hair follicular bulge SCs, suggesting the existence of a 
common cluster of epithelial SC genes. As we found, they also 
detected an elevated expression of the Sry gene in mouse 



limbal basal cells, associating it with increased proliferation. 
They proposed that it is involved in SC activation, maintaining 
proliferative capacity needed for expansion of precursor cell 
populations, and for wound healing [48]. Similarly, Figueira 
et al. [32] in a microarray analysis to identify phenotypic 
markers of human limbal SCs in fetal and adult corneas, 
detected that Wnt-4 was differentially overexpressed in fetal 
limbus compared with central cornea. Its expression was 
restricted to the basal and immediate parabasal limbal 
epithelium of both the adult and fetal corneas. They suggested 
that, since Wnt-4 functions in diverse developmental phases 
involved in common morphogenic events, it was not 
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surprising that this gene was expressed by the basal limbal 
epithelium that plays a crucial role in differentiation [32]. 
Wnt-4 overexpression, together with high levels of KRT15, 
KRT14, and P-cadherin in limbal basal epithelium cells, was 
in concordance with the molecular expression profile of 
stratified epithelial tissues. These data are in complete 
agreement with our RT 2 ~PCR-A results that confirm an 
upregulated expression for both Wnt and KRT15 molecules in 
limbal-derived epithelial cells. 

Analysis of our RT 2 -PCR-A data with IPA software also 
constructed 4 networks that are distinct from canonical 
pathways because they were generated de novo based on our 
input data. The resulting networks require further studies to 
find the most useful genes for defining a potential LESC 
profile. 

Conclusions — In conclusion, our study has led to the 
identification of novel molecules, CXCL12, ISL1, COL2A, 
NCAM1, ACAN, FOXA2, GJB1/Cnx32, and MSX1, that 
potentially could serve to recognize LESCs. Other markers, 
NCAM1 and GJB1/Cnx32 positively and PARD6A 
negatively, could be used to separate the stem cell-containing 
population of limbal epithelial cells derived from limbal niche 
cells grown in culture and intended for transplantation. 
Furthermore, the functional analysis of our results has 
provided a better understanding of the signaling molecular 
pathways associated with the progenitor-rich limbal 
epithelium. This knowledge potentially could give support to 
the design and development of innovative therapies with the 
potential to reverse corneal blindness arising from ocular 
surface failure due to LSCD. 

ACKNOWLEDGMENTS 

We thank the Barraquer Eye Bank of Barcelona for their 
support in providing human corneoscleral tissues. The authors 
thank Victoria Saez for her technical assistance and B. 
Bromberg (Certified Editor in Life Science of Xenofile 
Editing, www.xenofileediting.com) for his assistance in the 
final editing and preparation of this manuscript. This work 
was supported by Instituto de Salud Carlos III (CIBER-BBN 
CB06/0 1/003) and Junta de Castilla y Leon, Spain (Centro en 
Red de Medicina Regenerativa y Terapia Celular de Castilla 
y Leon; SAN673/VA/28/08). M. L-P. and S. G. were 
supported by fellowships from the Junta de Castilla y Leon. 
Presented in part as an abstract at the 80th Annual Meeting of 
the Association for Research in Vision and Ophthalmology 
(ARVO) May 2008, Fort Lauderdale, FL (IOVS 2008; ARVO 
E-abstract 6074 (0146-0404)), at the Stem Cells Europe 
Conference and Exhibition, September 2008, Amsterdam, 
The Netherlands; and at the TERMIS-NA 2008 Conference 
and Exposition, December 2008, San Diego, CA. 

REFERENCES 

1. Li W, Hayashida Y, Chen YT, Tseng SC. Niche regulation of 
corneal epithelial stem cells at the limbus. Cell Res 2007; 
17:26-36. [PMID: 17211449] 



2. Cotsarelis G, Cheng SZ, Dong G, Sun TT, Lavker RM. 

Existence of slow-cycling limbal epithelial basal cells that can 
be preferentially stimulated to proliferate: implications on 
epithelial stem cells. Cell 1989; 57:201-9. [PMID: 2702690] 

3. Schermer A, Galvin S, Sun TT. Differentiation-related 

expression of a major 64K corneal keratin in vivo and in 
culture suggests limbal location of corneal epithelial stem 
cells. J Cell Biol 1986; 103:49-62. [PMID: 2424919] 

4. Tseng SC. Regulation and clinical implications of corneal 

epithelial stem cells. Mol Biol Rep 1996; 23:47-58. [PMID: 
8983018] 

5. Dua HS, Saini JS. zuara-Blanco A, and Gupta P. Limbal stem 

cell deficiency: concept, aetiology, clinical presentation, 
diagnosis and management. Indian J Ophthalmol 2000; 
48:83-92. [PMID: 11116520] 

6. Cauchi PA, Ang GS. zuara-Blanco A, and Burr JM. A 

systematic literature review of surgical interventions for 
limbal stem cell deficiency in humans. Am J Ophthalmol 
2008; 146:251-9. [PMID: 18486098] 

7. Tsai RJ, Li LM, Chen JK. Reconstruction of damaged corneas 

by transplantation of autologous limbal epithelial cells. N 
Engl J Med 2000; 343:86-93. [PMID: 10891515] 

8. Koizumi N, Inatomi T, Suzuki T, Sotozono C, Kinoshita S. 

Cultivated corneal epithelial stem cell transplantation in 
ocular surface disorders. Ophthalmology 2001; 
108:1569-74. [PMID: 11535452] 

9. Kolli S, Ahmad S, Lako M, Figueiredo F. Successful clinical 

implementation of corneal epithelial stem cell therapy for 
treatment of unilateral limbal stem cell deficiency. Stem Cells 
2010; 28:597-610. [PMID: 20014040] 

1 0. Rama P, Matuska S, Paganoni G, Spinelli A, De LM, Pellegrini 

G. Limbal stem-cell therapy and long-term corneal 
regeneration. N Engl J Med 2010; 363:147-55. [PMID: 
20573916] 

1 1 . Shortt AJ, Seeker GA, Notara MD, Limb GA, Khaw PT, Tuft 

SJ, Daniels JT. Transplantation of ex vivo cultured limbal 
epithelial stem cells: a review of techniques and clinical 
results. Surv Ophthalmol 2007; 52:483-502. [PMID: 
17719371] 

12. Shortt AJ, Seeker GA, Rajan MS, Meligonis G, Dart JK, Tuft 

SJ, Daniels JT. Ex vivo expansion and transplantation of 
limbal epithelial stem cells. Ophthalmology 2008; 
115:1989-97. [PMID: 18554721] 

13. Chen Z, De Paiva CS, Luo L, Kretzer FL, Pflugfelder SC, Li 

DQ. Characterization of putative stem cell phenotype in 
human limbal epithelia. Stem Cells 2004; 22:355-66. [PMID: 
15153612] 

14. Lyngholm M, Vorum H, Nielsen K, Ostergaard M, Honore B, 

Ehlers N. Differences in the protein expression in limbal 
versus central human corneal epithelium-a search for stem 
cell markers. Exp Eye Res 2008; 87:96-105. [PMID: 
18571161] 

15. Schlotzer-Schrehardt U, Kruse FE. Identification and 

characterization of limbal stem cells. Exp Eye Res 2005; 
81:247-64. [PMID: 16051216] 

16. Takacs L, Toth E, Berta A, Vereb G. Stem cells of the adult 

cornea: from cytometric markers to therapeutic applications. 
Cytometry A 2009; 75:54-66. [PMID: 19051301] 



2115 



Molecular Vision 2011; 17:2102-2117 <http://www.molvis.org/molvis/vl7/a229> 

17. Pajoohesh-Ganji A, Stepp MA. In search of markers for the stem 

cells of the corneal epithelium. Biol Cell 2005; 97:265-76. 
[PMID: 15762848] 

18. Notara M, Daniels JT. Biological principals and clinical 

potentials of limbal epithelial stem cells. Cell Tissue Res 
2008; 331:135-43. [PMID: 17701219] 

19. Chen Z, Evans WH, Pflugfelder SC, Li DQ. Gap junction 

protein connexin 43 serves as a negative marker for a stem 
cell-containing population of human limbal epithelial cells. 
Stem Cells 2006; 24:1265-73. [PMID: 16424398] 

20. de Paiva CS, Chen Z, Corrales RM, Pflugfelder SC, Li DQ. 

ABCG2 transporter identifies a population of clonogenic 
human limbal epithelial cells. Stem Cells 2005; 23:63-73. 
[PMID: 15625123] 

21. Li DQ, Chen Z, Song XJ, De Paiva CS, Kim HS, Pflugfelder 

SC. Partial enrichment of a population of human limbal 
epithelial cells with putative stem cell properties based on 
collagen type IV adhesiveness. Exp Eye Res 2005; 
80:581-90. [PMID: 15781286] 

22. Qi H, Li DQ, Shine HD, Chen Z, Yoon KC, Jones DB, 

Pflugfelder SC. Nerve growth factor and its receptor TrkA 
serve as potential markers for human corneal epithelial 
progenitor cells. Exp Eye Res 2008; 86:34-40. [PMID: 
17980361] 

23. Qi H, Li DQ, Bian F, Chuang EY, Jones DB, Pflugfelder SC. 

Expression of glial cell-derived neurotrophic factor and its 
receptor in the stem-cell-containing human limbal epithelium. 
Br J Ophthalmol 2008; 92:1269-74. [PMID: 18723744] 

24. Corrales RM, Calonge M, Herreras JM, Saez V, Mayo A, 

Chaves FJ. Levels of mucin gene expression in normal human 
conjunctival epithelium in vivo. Curr Eye Res 2003; 
27:323-8. [PMID: 14562169] 

25. Corrales RM, Calonge M, Herreras JM, Saez V, Chaves FJ. 

Human epithelium from conjunctival impression cytology 
expresses MUC7 mucin gene. Cornea 2003; 22:665-71. 
[PMID: 14508262] 

26. Corrales RM, Galarreta DJ, Herreras JM, Calonge M, Chaves 

FJ. Normal human conjunctival epithelium expresses 
MUC13, MUC15, MUC16 and MUC17 mucin genes. Arch 
Soc Esp Oftalmol 2003; 78:375-81. [PMID: 12898407] 

27. Corrales RM, Stern ME, De Paiva CS, Welch J, Li DQ, 

Pflugfelder SC. Desiccating stress stimulates expression of 
matrix metalloproteinases by the corneal epithelium. Invest 
Ophthalmol Vis Sci 2006; 47:3293-302. [PMID: 16877394] 

28. Corrales RM, Villarreal A, Farley W, Stern ME, Li DQ, 

Pflugfelder SC. Strain-related cytokine profiles on the murine 
ocular surface in response to desiccating stress. Cornea 2007; 
26:579-84. [PMID: 17525655] 

29. De Paiva CS, Villarreal AL, Corrales RM, Rahman HT, Chang 

VY, Farley WJ, Stern ME, Niederkorn JY, Li DQ, Pflugfelder 
SC. Dry eye-induced conjunctival epithelial squamous 
metaplasia is modulated by interferon-gamma. Invest 
Ophthalmol Vis Sci 2007; 48:2553-60. [PMID: 17525184] 

30. Narayanan S, Corrales RM, Farley W, McDermott AM, 

Pflugfelder SC. Interleukin-1 receptor- 1 -deficient mice show 
attenuated production of ocular surface inflammatory 
cytokines in experimental dry eye. Cornea 2008; 27:811-7. 
[PMID: 18650668] 

31. Tong L, Corrales RM, Chen Z, Villarreal AL, De Paiva CS, 

Beuerman R, Li DQ, Pflugfelder SC. Expression and 



© 201 1 Molecular Vision 

regulation of cornified envelope proteins in human corneal 
epithelium. Invest Ophthalmol Vis Sci 2006; 47:1938-46. 
[PMID: 16639001] 

32. Figueira EC, Di GN, Coroneo MT, Wakefield D. The phenotype 

of limbal epithelial stem cells. Invest Ophthalmol Vis Sci 
2007; 48:144-56. [PMID: 17197527] 

33. Pellegrini G, Dellambra E, Golisano O, Martinelli E, Fantozzi 

I, Bondanza S, Ponzin D, McKeon F, De LM. p63 identifies 
keratinocyte stem cells. Proc Natl Acad Sci USA 2001; 
98:3156-61. [PMID: 11248048] 

34. Yang A, Schweitzer R, Sun D, Kaghad M, Walker N, Bronson 

RT, Tabin C, Sharpe A, Caput D, Crum C, McKeon F. p63 is 
essential for regenerative proliferation in limb, craniofacial 
and epithelial development. Nature 1999; 398:714-8. [PMID: 
10227294] 

35. Yoshida S, Shimmura S, Kawakita T, Miyashita H, Den S, 

Shimazaki J, Tsubota K. Cytokeratin 15 can be used to 
identify the limbal phenotype in normal and diseased ocular 
surfaces. Invest Ophthalmol Vis Sci 2006; 47:4780-6. 
[PMID: 17065488] 

36. Kasper M, Moll R, Stosiek P, Karsten U. Patterns of cytokeratin 

and vimentin expression in the human eye. Histochemistry 
1988; 89:369-77. [PMID: 2457569] 

37. Du Y, Chen J, Funderburgh JL, Zhu X, Li L. Functional 

reconstruction of rabbit corneal epithelium by human limbal 
cells cultured on amniotic membrane. Mol Vis 2003; 
9:635-43. [PMID: 14685149] 

38. Dua HS, Joseph A, Shanmuganathan VA, Jones RE. Stem cell 

differentiation and the effects of deficiency. Eye (Lond) 2003; 
17:877-85. [PMID: 14631392] 

39. Di Iorio E, Barbaro V, Ruzza A, Ponzin D, Pellegrini G, De LM. 

Isoforms of DeltaNp63 and the migration of ocular limbal 
cells in human corneal regeneration. Proc Natl Acad Sci USA 
2005; 102:9523-8. [PMID: 15983386] 

40. Adachi W, Ulanovsky H, Li Y, Norman B, Davis J, Piatigorsky 

J. Serial analysis of gene expression (SAGE) in the rat limbal 
and central corneal epithelium. Invest Ophthalmol Vis Sci 
2006; 47:3801-10. [PMID: 16936091] 

4 1 . Akinci MA, Turner H, Taveras M, Barash A, Wang Z, Reinach 

P, Wolosin JM. Molecular profiling of conjunctival epithelial 
side-population stem cells: atypical cell surface markers and 
sources of a slow-cycling phenotype. Invest Ophthalmol Vis 
Sci 2009; 50:4162-72. [PMID: 19324848] 

42. Akinci MA, Turner H, Taveras M, Wolosin JM. Differential 

gene expression in the pig limbal side population: 
implications for stem cell cycling, replication, and survival. 
Invest Ophthalmol Vis Sci 2009; 50:5630-8. [PMID: 
19608544] 

43 . Bian F, Liu W, Yoon KC, Lu R, Zhou N, Ma P, Pflugfelder SC, 

Li DQ. Molecular signatures and biological pathway profiles 
of human corneal epithelial progenitor cells. Int J Biochem 
Cell Biol 2010; 42:1142-53. [PMID: 20363357] 

44. Ding Z, Dong J, Liu J, Deng SX. Preferential gene expression 

in the limbus of the vervet monkey. Mol Vis 2008; 
14:2031-41. [PMID: 18989386] 

45. Jun AS, Liu SH, Koo EH, Do DV, Stark WJ, Gottsch JD. 

Microarray analysis of gene expression in human donor 
corneas. Arch Ophthalmol 2001; 119:1629-34. [PMID: 
11709013] 



2116 



Molecular Vision 2011; 17:2102-2117 <http://www.molvis.org/molvis/vl7/a229> 

46. Turner HC, Budak MT, Akinci MA, Wolosin JM. Comparative 

analysis of human conjunctival and corneal epithelial gene 
expression with oligonucleotide microarrays. Invest 
Ophthalmol Vis Sci 2007; 48:2050-61. [PMID: 17460260] 

47. Utheim TP, Raeder S, Olstad OK, Utheim OA, de La PM, 

Cheng R, Huynh TT, Messelt E, Roald B, Lyberg T. 
Comparison of the histology, gene expression profile, and 
phenotype of cultured human limbal epithelial cells from 
different limbal regions. Invest Ophthalmol Vis Sci 2009; 
50:5165-72. [PMID: 19578011] 

48. Zhou M, Li XM, Lavker RM. Transcriptional profiling of 

enriched populations of stem cells versus transient amplifying 
cells. A comparison of limbal and corneal epithelial basal 
cells. J Biol Chem 2006; 281:19600-9. [PMID: 16675456] 

49. Takacs L, Toth E, Losonczy G, Szanto A, Bahr-Ivacevic T, 

Benes V, Berta A, Vereb G. Differentially expressed genes 
associated with human limbal epithelial phenotypes: new 
molecules that potentially facilitate selection of stem cell- 
enriched populations. Invest Ophthalmol Vis Sci 2011; 
52:1252-60. [PMID: 21071743] 

50. Zlotnik A, Yoshie O. Chemokines: a new classification system 

and their role in immunity. Immunity 2000; 12:121-7. [PMID: 
10714678] 

51. Bourcier T, Berbar T, Paquet S, Rondeau N, Thomas F, 

Borderie V, Laroche L, Rostene W, Haour F, Lombet A. 
Characterization and functionality of CXCR4 chemokine 
receptor and SDF-1 in human corneal fibroblasts. Mol Vis 
2003; 9:96-102. [PMID: 12677165] 

52. Kitaori T, Ito H, Schwarz EM, Tsutsumi R, Yoshitomi H, Oishi 

S, Nakano M, Fujii N, Nagasawa T, Nakamura T. Stromal 
cell-derived factor 1/CXCR4 signaling is critical for the 
recruitment of mesenchymal stem cells to the fracture site 
during skeletal repair in a mouse model. Arthritis Rheum 
2009; 60:813-23. [PMID: 19248097] 

53. Zhuang Y, Chen X, Xu M, Zhang LY, Xiang F. Chemokine 

stromal cell-derived factor 1/CXCL12 increases homing of 
mesenchymal stem cells to injured myocardium and 
neovascularization following myocardial infarction. Chin 
Med J (Engl) 2009; 122:183-7. [PMID: 19187644] 

54. Ye J, Yao K, Kim JC. Mesenchymal stem cell transplantation 

in a rabbit corneal alkali burn model: engraftment and 
involvement in wound healing. Eye (Lond) 2006; 
20:482-90. [PMID: 15895027] 

55. Gao L, Joberty G, Macara IG. Assembly of epithelial tight 

junctions is negatively regulated by Par6. Curr Biol 2002; 
12:221-5. [PMID: 11839275] 

56. Ozdamar B, Bose R, Barrios-Rodiles M, Wang HR, Zhang Y, 

Wrana JL. Regulation of the polarity protein Par6 by TGFbeta 
receptors controls epithelial cell plasticity. Science 2005; 
307:1603-9. [PMID: 15761148] 

57. Thiery JP, Huang R. Linking epithelial-mesenchymal transition 

to the well-known polarity protein Par6. Dev Cell 2005; 
8:456-8. [PMID: 15809027] 



© 201 1 Molecular Vision 

58. Matic M, Petrov IN, Chen S, Wang C, Dimitrijevich SD, 

Wolosin JM. Stem cells of the corneal epithelium lack 
connexins and metabolite transfer capacity. Differentiation 
1997; 61:251-60. [PMID: 9203348] 

59. Wolosin JM. Cell markers and the side population phenotype 

in ocular surface epithelial stem cell characterization and 
isolation. Ocul Surf 2006; 4:10-23. [PMID: 16669522] 

60. Prochnow N, Dermietzel R. Connexons and cell adhesion: a 

romantic phase. Histochem Cell Biol 2008; 130:71-7. [PMID: 
18481075] 

61. Hirabayashi Y, Yoon BI, Tsuboi I, Huo Y, Kodama Y, Kanno 

J, Ott T, Trosko JE, Inoue T. Protective role of connexin 32 
in steady-state hematopoiesis, regeneration state, and 
leukemogenesis. Exp Biol Med (Maywood) 2007; 
232:700-12. [PMID: 17463168] 

62. Hirabayashi Y, Yoon BI, Tsuboi I, Huo Y, Kodama Y, Kanno 

J, Ott T, Trosko JE, Inoue T. Membrane channel connexin 32 
maintains Lin(-)/c-kit(+) hematopoietic progenitor cell 
compartment: analysis of the cell cycle. J Membr Biol 2007; 
217:105-13. [PMID: 17629738] 

63. Crossin KL, Krushel LA. Cellular signaling by neural cell 

adhesion molecules of the immunoglobulin superfamily. Dev 
Dyn2000; 218:260-79. [PMID: 10842356] 

64. Watanabe M, Kobayashi H, Rutishauser U, Katar M, Alcala J, 

Maisel H. NCAM in the differentiation of embryonic lens 
tissue. Dev Biol 1989; 135:414-23. [PMID: 2776975] 

65. Filiz S, Dalcik H, Yardimoglu M, Gonca S, Ceylan S. 

Localization of neural cell adhesion molecule (N-CAM) 
immunoreactivity in adult rat tissues. Biotech Histochem 
2002; 77:127-35. [PMID: 12229933] 

66. Dua HS, Gomes JA, Jindal VK, Appa SN, Schwarting R, Eagle 

RC Jr, Donoso LA, Laibson PR. Mucosa specific 
lymphocytes in the human conjunctiva, corneoscleral limbus 
and lacrimal gland. Curr Eye Res 1994; 13:87-93. [PMID: 
7908866] 

67. Vantrappen L, Geboes K, Missotten L, Maudgal PC, Desmet V. 

Lymphocytes and Langerhans cells in the normal human 
cornea. Invest Ophthalmol Vis Sci 1985; 26:220-5. [PMID: 
3882607] 

68. Aandahl EM, Torgersen KM, Tasken K. CD8+ regulatory T 

cells-A distinct T-cell lineage or a transient T-cell phenotype? 
Hum Immunol 2008; 69:696-9. [PMID: 18817824] 

69. Clevers H. Wnt/beta-catenin signaling in development and 

disease. Cell 2006; 127:469-80. [PMID: 17081971] 

70. Lefebvre V, Dumitriu B, Penzo-Mendez A, Han Y, Pallavi B. 

Control of cell fate and differentiation by Sry-related high- 
mobility-group box (Sox) transcription factors. Int J Biochem 
Cell Biol 2007; 39:2195-214. [PMID: 17625949] 

7 1 . Lowry WE, Blanpain C, Nowak JA, Guasch G, Lewis L, Fuchs 

E. Defining the impact of beta-catenin/Tcf transactivation on 
epithelial stem cells. Genes Dev 2005; 19:1596-61 1. [PMID: 
15961525] 



Articles are provided courtesy of Emory University and the Zhongshan Ophthalmic Center, Sun Yat-sen University, P.R. China. 
The print version of this article was created on 10 August 2011. This reflects all typographical corrections and errata to the 
article through that date. Details of any changes may be found in the online version of the article. 

2117 



